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In the variational approach to quantum statistics, a smearing formula 
describes efficiently the consequences of quantum fluctuations upon an inter- 
action potential. The result is an effective classical potential from which the 
QQ ', partition function can be obtained by a simple integral. In this work, the 

0^ I smearing formula is extended to higher orders in the variational perturbation 

' theory. An application to the singular Coulomb potential exhibits the same 

^ ■ fast convergence with increasing orders that has been observed in previous 

^ . variational perturbation expansions of the anharmonic oscillator with quartic 

, potential. 

in . 

^ ■ I. INTRODUCTION 

^ : 

Q ■ The variational approach to quantum statistics, initially based on the Jensen-Peierls 

^ ! inequality for imaginary-time path integrals |1|J^, yielded crude upper bounds for the free 

^ I energy of many quantum mechanical systems at all temperatures and coupling strengths. 

OS I By abandoning the inequality, the approach has been extended to a systematic, fast con- 

^ ' vergent variational perturbation theory within which the original approach is just a 
first-order approximation. Thermodynamic and some local quantities can now be evaluated 

! to any desired accuracy, starting out from an ordinary perturbation expansion of arbitrary 

! order. 

^ : 

CT' A particularly attractive feature of the original variational approach was the existence 

'i> ■ of a smearing formula in the form of a Gaussian convolution integral which compactly 

^ , accounts for the effect of quantum fluctuations upon the interaction potential and other 

^ I local quantities This formula was applicable to some classes of singular potentials 



such as the Coulomb potential There is a definite need for such a formula in higher 

orders of variational perturbation theory, which so far was based on Feynman diagrams, 
thus being limited to polynomial interactions. The purpose of this paper is to derive the 
desired higher-order smearing formula. This will be done in Section III after a brief review 
of variational perturbation theory in Section II. An application to the Coulomb potential 
is given in Section IV, where the effective classical potential is calculated to second order 
in the variational perturbation theory. Its zero-temperature limit yields in Section V a 
variational perturbation expansion for the ground state energy up to second order. Section 
VI reproduces this result by a direct variational treatment of the Rayleigh-Schrodinger 
perturbation expansion, and carries it to third order to demonstrate the fast convergence of 
the variational perturbation expansion. 
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II. REVIEW OF VARIATIONAL PERTURBATION THEORY 



Consider a quantum mechanical point particle of mass M moving in a one-dimensional 
time-independent potential V{x). Its thermodynamic partition function is given by the 
imaginary-time path integral 

Z = ^Pa;(r)exp|-i^[x(r)]| , (1) 



with the euclidean action 

hi3 



A[xir)] = jdr[fxirr + Vixir))\ , (2) 



and the abbreviation (3 = l/fc^T. The paths x{t) satisfy the periodic boundary condition 
a;(0) = x{h(3). Following Feynman we decompose the path integral for the partition 
function ([^) into an ordinary integral over the time-averaged position 

If 

Xq = x = — / c/rx(r) , (3) 

and a remaining path integral over the fluctuations 

5x{t) = x(r) — xq (4) 
around Xq. Thus we rewrite (||) as an integral 

+ 00 



dxo (5) 

over a local partition function Z^" which is defined by the restricted path integral 

= ^Px(r)5(x-Xo)exp|-i^[x(r)]| , (6) 

with the notation 



The free energy associated with the local partition function (||) is defined as the effective 
classical potential |^ 

F<=ff.=i(xo) = -ilogZ-° (8) 

which accounts for the effects of all quantum fluctuations. 

In order to calculate V'^^''^^{xo), we decompose the euclidean action (0) into a sum 
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A\x{T)]=A^Mr)]+^^A<r)], (9) 

where the first term is the action of a harmonic oscillator centered around xq with an 
undetermined local trial frequency fi(xo), 

ft/3 

A^^W)] = I dr^^xirY + fn^xo) [x{t) - Xof] , (10) 



and the second term is the remaining interaction 

m<r)]= jdrV,ll{x{r)) (11) 



of the potential difference 

V.2{x)=V{x)-^n\x,){x-x,f. (12) 

With this decomposition, the local partition function (|^) may be expanded in powers of the 
interaction potential (|12]) around the local harmonic partition function 



Z^- = jvx{r)~5{x-xo)e^Y>[-\A^Mr)]] • (13) 

The expansion reads 

{hp np hf3 

1 - ^ / dn{v^^{x{nW^^ ^wl'^'^'l ^^2(^nn^(ri))v;:n^(^2)))^° 


ft/3 ft/3 ft/3 1 

'Wl^^^'l I ^^3(Knn^(ri))KS'(^(r2))l^iS'(^(r3)))?^" + . . . , (14) 
J 

where the xo-dependent expectation values (Fi(x(ri)) . . . F„(x(r„)))^° are correlation func- 
tions of the local harmonic trial system: 

(Fi(x(ri))---F„(x(r„)))?^°- ^ 



X ^Px(r)Fi(x(ri))---F„(x(r„))5(x-a;o)exp|-i^^°[x(r)]| . (15) 

The correlation functions can be decomposed into connected ones via the standard cumulant 
expansion yielding for the effective classical potential V^^'^^^Xq) the following pertur- 
bation expansion 

ft/3 ft/3 ft/3 

V^^^'^'ixo) = F^^ + ^/rfri(l^:n^(ri)))?i°-^/rfn/rfr2(l^:n^(n))V;-'(^(r2)))o°c 

'^00 

ft/3 ft/3 ft/3 

+ ^ f dn I dr, I c/r3(V;Sl'(x(ri))K::nx(r2))V;Sn^(r3)))?^!. + • • • • (16) 
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The first term on the right-hand side is the free energy of the local harmonic partition 
function 

The second term contains the local harmonic expectation value of the potential for which 
there exists the above mentioned smearing formula which we want to extend in this work. 
The cumulant in the third term is given by the following combination of expectation values: 

{V-2{<n))V.2{x{T,)))l^^^ = {VrMri))Vrj{x{T,W^> - {V^ti^inW^^ {VrMr2)m' ■ (18) 

By construction, the effective classical potential V'^^''^^{xo) in (|1^) does not depend on 
the choice of the frequency fl{xo) in the trial action (|10]). However, when truncating the 
infinite sum (|TBp after the iVth order, we obtain an approximation W^{xq) for the effec- 
tive classical potential V'^^''^^{xq) with an f2(xo)-dependence, which decreases with increasing 
order N of the expansion. With the expectation that the optimal truncated sum W^{xo) de- 
pends minimally on the frequency Q{xo), we therefore determine fi(xo) from the extremality 
condition 

= (19) 

If this has no solution, we demand as the next-best condition of minimal dependence on 
fi(xo) iii 

^^^"(^o) - (20) 

The result is called the optimal frequency Qn{xo) of order A^. It yields the truncated sum 
Wn{xo) = W^'^^^°\xq) which represents the desired A^th order approximation to the ef- 
fective classical potential V^^'^^^xq). The first-order approximation Wi{xo) coincides with 
the original variational result of Feynman and Kleinert which satifies the Jensen- Peierls 
inequality and guarantees the existence of an extremum (|T^). 

The accuracy of the approximate effective classical potential W^ixo) can be assessed 
by the following considerations In the limit of high temperatures where (3—^0, the 
approximation Wn{xo) becomes exact for any A^: 

lmiiy;v(xo) = V^(xo). (21) 

At low temperatures, where /3 oo, we obtain from (^ and (|^) an approximation to the 
free energy in form of an integral over the time-averaged position xq 




+ 00 

dxoexp[-(3WN{xo)]} , (22) 



whose integrand is centered sharply around the minimum x™™ of Wi\f{xo). Performing this 
integral in the saddle point approximation yields an A^-th order approximation E^^ for the 
ground state energy E^^^ of the quantum system. 
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E^^> = min lim WNixo) . (23) 

The same approximation to the ground state energy can also be obtained by a variational 
resummation of the Rayleigh-Schrodinger perturbation series for E'^^\ This will be shown 
in Section V for the ground state energy of the Coulomb potential up to the order N = 2. 



III. EVALUATION OF PATH INTEGRALS 

In order to calculate the different terms in the variational perturbation expansion (f6|), 
we must find efficient formulas for evaluating local correlation functions of the type (|15|). 
For this we observe that, by fixing of the temporal average ai x = xq in the path integral 



in (16), the zero Matsubara frequency = is removed from the Fourier decomposition of 



the periodic paths 

oo 

<r) = a;o + E (a^-e'"'"" + ^*m(^~''""'") , = 27im/hp . (24) 



m=l 



In fact, the restricted integration measure §Vx{t)S(x — xq) in (^) may be decomposed 
into a product of ordinary integrals over real and imaginary parts and of the Fourier 
components Xm according to 1^: 



ivx{r)~5{x~x,) =f[{ f dxZ f rfxr^^] . (25) 

m=l ^ J 

The zero-frequency component Xq remains unintegrated. With this decomposition, the basic 
local pair correlation function of the fluctuations 6x{t) in (^) can immediately be calculated 
from ([l5|) as a Matsubara sum without the zero mode: 

GS(ry) . iHrMS - £ (26) 

Performing the sum yields the explicit result 

r-ofr r') = ^ / cosh[n{xo)\r-r'\-hpn{xo)/2] _ 2 \ 

^ ^ ' ^ 2Mn{xo) \ smh[hi3n{xo)/2] h/3n{xo) J ■ ^ ^ 

The first term is the ordinary oscillator correlation function of frequency Q 

^^J^^;^. (28) 

while the last term subtracts the zero mode which is absent in (^). This absence has the 
important consequence that 

|rfrG?^»(r,r')=0. (29) 
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Using (^), the expectation values in (|T6D can easily be calculated for a polynomial 
potential using Wick's contraction rules, by which the expectation values can be reduced 
to sums over products of pair correlation functions (r, r') . In order to abbreviate the 
notation and to emphasize the dimension (length)^ of these quantities, we shall denotes the 
local Green functions in ( pTf ) from now on by a^^/(a;o). The harmonic expectation value of 
any odd power n in the fluctuation variable 6x{t) is zero. For even n, the Wick expansion 
reads 



Xo 



I n P 



(30) 



\fc=i 



where the sum runs over all (n — 1)!! pair contractions. For an exponential. Wick's rule 
implies 



exp 





xo 


i J dTj{T)5x{T) 


^ = exp 





n 







(31) 



In the special case j^r) = I]fc=itifc5(T — r^), we obtain the important formula for the ex- 
pecation value of a product of exponentials 



\k=l 



xo 



^iUkSx{Tk) 



exp 



n n 

'9 E E Ukal^^^Xxo)Uk' 



k=l k'=l 



(32) 



After Fourier-decomposing the functions Fi(x), . . . , F„(x) in (p^), formula ( |32D yields di- 
rectly the desired smearing formula [l^ 



(Fi(x(ri))---F„(a;(r„))) 



xo 



Yl / dxkFk{xk) 



X 



(27r)«Det a^^^X^o) 



exp 



n n 

- 9 E E ^^fc «r,r,, M 5Xk> 



k=l k' = l 



(33) 



where a~^^, (a;o) denotes the inverse of the sjTumetric n x n-matrix a^^^^, (xq). This smearing 
formula determines the different harmonic expectation values in the variational perturbation 
expansion (|16D as convolutions with Gaussian functions. 

For n = 1, the smearing formula (|33|) reduces to the previous one 



dxiFi(xi) 



1 



27ra^(a;o) 



: exp 



2a2(xo) 



(34) 



where a^(xo) denotes the r-independent diagonal matrix element a^^(xo). For polynomials 
Fi{x), the smearing formula (|33| ) reproduces Wick's rule: Odd powers in 6x{t) have van- 
ishing local correlation functions, whereas even powers result in (^), which for coinciding 
times Tk reduces to 



([Mr.)re = (^-1)!! a"(xo). 



(35) 
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For two functions Fi{x) and F2{x), our smearing formula (|33|) reads, more explicitly, 

1 



hoo +00 



{Fi(x(ri))F2(x(r2)))^° = J dxi J dx2 Fi{x,)F2{x2) 



—00 —00 



X exp 



a^(xo)(xi - xoY - 2al^^,^{xo)ixi - xo){x2 - xq) + a^{xo){x2 - xqY 



2[a^ixo) - ai^^^ixo)] 
Specializing ^2(0;) to the square of the function 6x, we obtain the useful rule 



(36) 



Fi(x(ri)) MT2)fy^ = a\xo) 



a4(xo) 



XQ 

n 



+ 



(37) 



which reduces the smearing procedure for different times ti and T2 to corresponding ones at 
equal times Ti = T2. With this we immediately yield, for instance. 



[6x{n)f[6x{T2)f)^ = a\x,)+2at^^,^{xo). 



(38) 



In three dimensions, the trial potential contains a 3 x 3-frequency matrix Qij depending 
on the time-averaged position xq and reads 



M ^ 



(39) 



while the interaction potential (|1^) becomes 



M 



(40) 



Assuming the potential to depend only on r = |x|, i.e. ^(x) = v{r), the frequency matrix 
possesses only two invariant matrix elements, a longitudinal one f2i(ro) and a transversal 
one nr(ro) The interaction potential (^) can then be decomposed into a longitudinal 

and a transversal part according to 



(41) 



so that (HO) may be rewritten as 



M 



{nliro) Ml + nUro) Ml} , 



(42) 



with obvious definitions of the longitudinal and transverse projections ^x^, and (5xt of the 
fluctuations (5x. 

To first order, the anisotropic generalization of the smearing formula (|34D reads |^J^ 
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exp 



2al 



2al 



(43) 



For the squares of transverse and longitudinal fluctuations, this generalizes (|35D with n = 2 
to 



2a 



T ! 



The second-order smearing formula (p6D becomes in three dimensions 

+ 00 +00 



X 



exp 



X exp 



'' al{xiL - ro)^ - 2al^^^^{xiL - ro)(x2L - tq) + a|(x2L - ro)^ 1 



so that the rule { ^7\j for expectation values is generalized to 



2ai 



a 



1 - 



''TT1T2 



(i^i(x(ri))) 



ro 



TtiT2 



Fi(x(n))[5x(ri)]T>j;'^,^^ 



a 



LTir2 



(i^i(x(ri)))- 
, ,Fi(x(ri)[5x(ri)]L 



''LT1T2 



(44) 



(45) 



(46) 



(47) 



Specializing -Fi(x) to a quadratic function, we obtain the corresponding generalizations of 



[5x(ri)]^ Mr2)]iyi^^^ 
[<5x(ri)]?^[5x(r2)]i);;^^^ 
[<5x(rO]i [5x(r2)]i);;'^^^^ 



4a^ + 4a^ 



Tt\T2 ' 



Lt\T2 ' 



(48) 
(49) 
(50) 



IV. APPLICATION TO COULOMB POTENTIAL 

Let us demonstrate the use of the new smearing formulas by calculating the effective 
classical potential of the three-dimensional Coulomb potential 

v{r) = (51) 
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up to the second order in the variational perturbation expansion, thus going beyond the 
known first-order results in 0,0. For harmonic and Coulomb potentials we express the 
Coulomb potential (^) as a "proper-time" integral 



1 

- = An 
r 



(27r)3 



dae 



-crk'^— ikx 



(52) 



where a has the dimension (length)^, and find the expectation value 



I2al 



TT 



dX 



al - ai)A2 + al 



exp 



2al 



(53) 



From a straight-forward three-dimensional extension of (|TB|), the first-order variational ap- 
proximation to the effective classical potential is then 0,^] 



2 smh[n(3nT/2] , 

^'^^ = p^'^ npnr/2 +^^"- npnj2 



1 . smh[ni3nL/2] M 



'-^{2nWT + ^>l} 



/2ai 

TT 



dX' 



1 



exp 



2ai 



(54) 



where we have ommitted the argument tq from all functions on the right-hand side. A 
similar expression was derived in the isotropic approximation fli = with the help of 
Gaussian wave packets in the context of plasma physics for the purpose of faster molecular 
dynamics simulations [p!l| -[T3[|. 

]) into the rules 



By inserting the "proper-time" integral for the Coulomb potential 
\) and (0), we find 



— -— [5x{t2)]t ) 



ro 



'2a| 

71 



X 



dX exp 



2a| 



2al 



'2a| 

71 



al-al)X^ + al [{a^ - al)X^ + al]" 
1 . 



dX exp 



2al 



(55) 
(56) 



Note that these results are also special cases of the general expectation value 



1 v 



X exp 



27r2 



da 



exp 1 






2 







d^xF(x) 



2a)x: 



[ai + 2a){xL - ro)2 + 2al^^^^ro{xL - ro) 



2[a^ 



Ttit2 



+ 2a^crl 



2K 



+ 2a?crl 



(57) 



which gives us the additional local expectation value 
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X Ti X Ts 



ro 



2 

TT 



dai / (icT; 



X 



[al + 2(Ti][al + 2a2 



exp 



[a| + 2ai][a| + 2a2] - 



''Lrir2 



[ai + 2ai][ai + 2a2] 



(58) 



From these smearing results we calculate the connected correlation functions of the interac- 
tion potential (^Op according to the cumulant law (|TB|), and insert these into (|T^ to obtain 
the second-order effective classical potential 



TX 



M2[21]|4 + n|/|] 



Ah 



at[{a'^-al)X' + al]l 



exp 



2al' 



|x(ri)| |x(r2) 



with the abbreviation 



4 + h'^p'^n^L - 4 cosh hpnT,L + ^P^t,l sinh hpnT,L 



SpM^n^ L sinh^ [h[3nT,L/2] 



(60) 



which is a function of tq of dimension (length)^. After extremizing (E3) and (pUI) with re- 



spect to the trial frequencies VLt^^l according to ([T9|), which has to be done numerically 
we obtain the first- and second-order approximations for the effective classical potential of 
the Coulomb system. The isotropic approximations VLt = VLl m the first and second order 
are plotted in Fig. 1 for various temperatures. The second-order curves lie all below the 
first-order ones, and the difference between the two decreases with increasing temperature 
and increasing distance from the origin. Figure 2 shows exemplarily that the anisotropic 
approximation slightly deviates from the isotropic one. The difference between both is only 
visible for intermediate distances from the origin. 



V. ZERO-TEMPERATURE LIMIT 



In order to check our results we take (|54D and (|59| ) to the limit T ^ 0, where W^{rQ) 
and 1^2^ (ro) reduce at rg = 0, according to (p3D , to the ground state energy of the Coulomb 
system. At ro = 0, the frequencies are isotropic Vt^ = = Q for symmetry reasons, thus 
simplifying the evaluation (Q) and (^91) . Taking into account the low temperature limit of 
the two-point correlations (P^) 

^'^''r.r.A^o) 2Mfi(Xo) ' ^ ^ 

we immediately deduce for the first order approximation (|5^) at xq = with Q = f2(0) the 
limit 
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-0.2 



-0.4 



-0.6 



-0.8 



V{ro) 



0.5 



1.5 



2 



2.5 



3.5 



FIG. 1. Isotropic approximations to the effective classical potential of the Coulomb system in 
the first (lines) and second order (dots). The temperatures are 0, 0.1, 0.01 and cxd from the top 
to the bottom in atomic units. The high temperature limit is the same for all approximations 
WN{ro). 



Ef\n) = lim <(0) = hn - -^J^e' 



(62) 



To second order, the limit is more involved. Performing the integrals over ai and (J2 in 
l), we obtain with (IIBD the connected correlation function 



1 



X Ti X Ts 



(0) vra^ (0) 



arctan 



anr2(0) 



Txa 



T1T2 



:o)- 



(63) 



Inserting here the zero-temperature limit (|6T1) , we can integrate these expressions over the 
imaginary times ti,T2 G [0, hP], and find for large /3 



ft/3 hl3 



n,c 



AM 



2 

IT 



X 



e^^^ arcsin ^ 1 - e-^^^^ + ^ log a{f3) - I [log a{(3)f - ^ [ du^ 

2 8 2 J L 

a{l3) 



logM 



+ U 



(64) 
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FIG. 2. Isotropic and anisotropic approximations to the effective classical potential of the 
Coulomb system in the first and second order at the temperature 0.1 in atomic units. The lowest 
line represents the high temperature limit in which all isotropic and anisotropic approximations 
coincide. 



with the abbreviation 



1 _ _ Q~2n/3n 
- 1 + VI - e-2?i?^ ■ ^^^^ 
Inserting this into (|59|) and going to the hmit /3 — oo, we find the ground state energy 



E^m = Ito »?(0) = j-^na - ^^—e' - - (^IH- log2 - -J ^e«. (66) 

Postponing for a moment the extremization of (|62D and (^) with respect to the trial 
frequency f2, let us first rederive this result from a variational treatment of the ordinary 
Rayleigh-Schrodinger perturbation expansion for the ground-state energy. 



VI. VARIATIONAL TREATMENT OF RAYLEIGH-SCHRODINGER 
PERTURBATION EXPANSION 

According to the rules layed out in 0, we must first calculate the ground-state energy 
for the Coulomb potential in the presence of a harmonic potential of frequency uj\ 
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Kux(x) = y^V-^. (67) 

After this, we make the trivial replacement uj v^T^ + oJ^— TF, and reexpand the energy 
in powers of the difference uo"^ — D?, considering this quantity as being of the order and 
truncating the reexpansion accordingly. At the end we set = 0, since the original Coulomb 
system contains no oscillator potential. Note that this limit is equivalent to a strong-coupling 
limit of (p7D with respect to e^. The result of this treatment will coincide precisely with the 
expansions (|^) and (|66D , respectively. 

The Rayleigh-Schrodinger perturbation expansion of the ground state energy E^^{(jj) for 
the potential (^) requires the knowledge of the matrix elements of the Coulomb potential 
(pT|) with respect to the eigenfunctions ipn,i,m{f','&,y^) of the harmonic oscillator with the 
frequency uj 0: 

27V n cso 2 

Vn,i,m;n',i',m' = j dif j smi^ J dr tp^^i „^{r, 1^ , (f) -^^n' ,1' ,m' {r, 1^ , (f) , (68) 





2n\ 



T{n + l + 3/2) 



(69) 



Here n denotes the radial quantum number, L'^{x) the Laguerre polynomials [14 , and 
Yi^jjii'^,'^) the spherical harmonics obeying the orthonormality relation 



2n TV 

d^ f d^sm^Y;^{^,if)Y,,^,{^,<f) = 5^,5„ 



(70) 







Inserting ( |5UD into (pSD, and evaluating the integrals with Eq. (2.19.14.15) in Ref. |T3 
find 



we 



2 MujT{l + l)T{n + l/2) T{n' + 1 + 3/2) 
~^ r(/ + 3/2) \ n\n'\r{n + 1 + 3/2) 

X3F2 (^-n', / + 1, / + ^, i - n; 1^ 5;,// 6m,m' , 



(71) 



with the generalized hypergeometric series |14] 



3^2(^1, aa, as; Pi, P2; x) = J2 



(ai)fc(a2)fc(a3)fc x'' 



(72) 



and the Pochhammer symbol {ot)k = r(a + k)/T{a). 

These matrix elements are now inserted into the Rayleigh-Schrodinger perturbation ex- 
pansion for the ground state energy 
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jp / \ _ T7 I I K),0,0;n,Z,mKi,i,m;0,0,0 ^0, 0,0;n, I, mVn, I, m;0, 0,0 

^a.ux[^^) — ^0,0,0 + »^0,0,0;0,0,0 + 2^ ^ 2^ V^0,0,0;0,0,0— j— =; 72~ 

n,l,m ^0,0,0 ^n,l,m [-i^0,0,0 " ^n,l,m\ 

_j_ ^ ^ yO,0,0;n,l,myn,l,m;n', I', m'Vn', I', m' -,0,0,0 _|_ ^yg-^ 
n,l,mn',l',m' ["^0,0,0 ~ -Sn,Z,m] [-^0,0,0 ~ -Sn',Z',m'] 

the denominators containing the energy eigenvalues of the harmonic oscillator 

En,i,m = hu;(2n + l + ^^ . (74) 

The primed summations in run over all values of the quantum numbers n, I = 

—oo, . . . , +00 and m = — /, . . . , +/, excluding those for which the denominators vanish. 
For the first three orders we obtain from ([7I|)-([7^ 



, , 3^ 2 Muj n 4 / , 7r\ M , . , , 

= _S„ ^ _ _ (1 + ,og 2 - -) ^e' - c + . . . (76) 



with the constant 

1 f ^ 1 ■ 3 ■ . . . ■ (2n - 1) 1 ~ ^ 1 ■ 3 ■ . . . ■ (2n - 1) 



c 



7r3/2\J^ 2-4-...-2n n^n + 1/2) J.^^i 2 ■ 4 ■ . . . ■ 2n 



1 ■ 3 ■ . . . ■ {2n' - 1) 3F2 {-n', / + 1, |; / + |, | - n; 1 
2 ■ 4 ■ . . . ■ 2n' nn' (n + 1/2) 



0.0318. (76) 



The variational reexpansion procedure described after Eq. ( p7D replaces in the first term by 
fi(l — 1)^/^, to be expanded in the second 1 up to third order as 1 — 1/2 — 1/8 — 1/16 = 5/16. 
Correspondingly the term 3a;/2 becomes 15/32a;. The factor cj^/^ goes over into fi^/^(l — 
1)^/^, to be expanded to second order in the second 1, yielding f2^/'^(l — 1/4 — 3/32) = 21/32. 
The next term in (|75| ) happens to be independent of u and needs no reexpansion, whereas 
the last term remains unchanged since it is of highest order. In this way we obtain from 
([75|) the third-order variational expression 



We are now ready to optimize successively the expansions of first, second, and third order 
(|6^), (|6^) and ( ffTl) with respect to the trial frequency Q. From the extrema we find the 
frequencies 



1^1 = 1^2 = -^, ^3 = C ^ (78) 

where c' ~ 0.7254 is the largest of the three solutions of the cubic equation 15c'^ — 21c'^/7r + 
16c. The corresponding approximations to the ground state energy are 

EP{n^) = -7N^, (79) 
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with the constants 



7i 



— ^ 0.424, 

3n 



72 



5 + 41og2 



0.474 , 



TT 



73 ^ 0.490, 



which quickly approaching the exact value 7ex = 0.5, as shown in Fig. 3. 



1 



3 N 



In 



-0.42 
-0.44 
-0.46 
-0.48 
-0.5 




3 N 



exact 



(80) 



FIG. 3. Rapid approach of the variational approximations to the ground state energy E^\Qn) 
to the correct ground state energy —0.5 (in natural energy units Me'^/h^). The right-hand loga- 
rithmic plot shows a slope —1. 



VII. SUMMARY AND OUTLOOK 

In this work we have extended the rules for calculating higher orders in variational 
perturbation theory from polynomial to nonpolynomial interactions. The effective classical 
potential of a quantum mechanical system is obtained from an extension of the known 
first-order smearing formula, and involves certain convolutions with Gaussian functions. As 
an example, we have applied the higher-order smearing formula to the Coulomb system. 
We have illustrated the fast convergence of the variational perturbation expansion even 
for such a singular potential. The new smearing formula will help improving the existing 
first-order variational results for partition function and density matrix, also in dissipative 
quantum systems |T^|T^. It will also be of use in treating field theories with nonpolynomial 
interactions such as Sine-Gordon and Liouville theories. 
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